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Dislocation-induced spin tunneling in Mn 12 acetate 
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Comprehensive theory of quantum spin relaxation in Mni2 acetate crystals is developed, that takes 
into account imperfections of the crystal structure and is based upon the generalization of the Landau- 
Zener effect for incoherent tunneling from excited energy levels. It is shown that linear dislocations 
at plausible concentrations provide the transverse anisotropy which is the main source of tunneling 
in Mni2. Local rotations of the easy axis due to dislocations result in a transverse magnetic field 
generated by the field applied along the c-axis of the crystal, which explains the presence of odd 
tunneling resonances. Long-range deformations due to dislocations produce a broad distribution 
of tunnel splittings. The theory predicts that at subkelvin temperatures the relaxation curves for 
different tunneling resonances can be scaled onto a single master curve. The magnetic relaxation in 
the thermally activated regime follows the stretched-exponential law with the exponent depending 
on the field, temperature, and concentration of defects. 

PACS numbers: 75.45. +j, 75.50.Tt 



I. INTRODUCTION 

Physical properties of M1112 acetate crystals have been 
subject of intensive investigation in the last years. The 
crystals have been chemically synthesized by LisEl in 1980, 
who found that they had a centered tetragonal struc- 
ture with *iri= 17.319A and c = 12.388A as lattice 
parametersira At the beginning of 1990's Sessoli et al 
of the Florence group of chemists (see, e.g., Refs. ||, H, 
||) established that Mni 2 clusters located at the sites of 
the crystal lattice have spin 10. Sessoli et al also demon- 
strated that Mni2 crystals are characterized by a record 
high uniaxial magnetic anisotropy along the c-direction, 
that makes the 65K energy barrier between the spin- 
up and spin-down states. The explosion of the inter- 
est of, physicists to Mni2 acetate came after Friedman 
et alu measured its spectacular stepwise magnetic hys- 
teresis and explained it by resonant spin tunneling (see 
also the followup experiments, Refs. [?], ||). Numerous 
experiments on Mni2 performed since 1996 uncovered a 
number, of other interesting phenomena. s.uch as memory 
effects,EI non-exponential relaxation,Lnl30 and a peculiax 
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crossover between tjifiana, 
predicted by theory. 



Uniaxial spin Hamiltonian, in the field parallel to the 
anisotropy axis, H = ~DS Z — H Z S Z , has pairs of degen- 
erate levels for 



H z = kD, k = 0,±1,±2,...,±(25- 1) 



(1) 



For a sample initially magnetized in the negative z- 
direction, Mni2 molecules occupy spin states with a neg- 
ative magnetic quantum number m. When the magnetic 
field is applied in the positive z-direction, the molecules 
eventually relax to the states with positive m. It may 



occur via thermal activation over the anisotropy barrier 
or through quantum tunneling between the states on dif- 
ferent sides of the barrier (see Fig. |l|). The latter process 
adds to the thermal activation when the levels on the two 
sides of the barrier are on resonance, that is, at H z = kD. 
Consequently, at H z — kD the magnetic relaxation of the 
crystal is faster than for the off-resonance values of the 
neld.Q 

There exist two macroscopic experimental approaches 
to the study of spin tunneling in Mni2 (see, e.g., Refs. 

1 0, 1 n l§, 6 i§ 0i> U SUSS and 

references therein). Both clearly demonstrate the effect 
of resonant spin tunneling. In the first approach, one 
measures the magnetic relaxation in the crystal at fixed 
temperature and magnetic field. The theory of such a 
relaxation works out tunneling transitions between lev- 
els on two sides of the energy barrier and spin-phonon 
transitions between the levels on one side of the barrier. 
These processes are described by the density-matrix for- 
malism suggested in Ref. |l3| and further developed in 
Refs. [33], |(| |34], |35|. In the second approach, one sweeps 
magnetic field through the resonant value and measures 
the fraction of molecules that change their magnetic mo- 
ments. For the ground-state tunneling, the theory of 
such process has been., djazekiped-along the lines of the 
Landau-Zener effect .EjO'ESOO'EI A nice feature of the 
field-sweep approach is that the fraction of molecules that 
change the direction of the spin depends on the tunneling 
level splitting and on the sweep rate but is insensitive to 
the dissipation in the experimental limit of small split- 
ting. 

While the uniaxial Hamiltonian allows to compute res- 
onant fields from independent measurements of D, it does 
not explain why spin tunneling in Mni2 actually occurs. 
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FIG. 1. Spin energy levels of a Mni2 Ac molecule for 
H x = and H z = D corresponding to the first resonance, 
k = 1. 



This is because TC = —DS Z — H Z S Z conserves 5*2 and 
thus does not allow any transitions between different m. 
To obtain tunneling, one should include in the Hamil- 
tonian the terms which do not commute with S z . The 
tetragonal symmetry of the crystal does not allow trans- 
verse anisotropy terms which are quadratic on the spin 
operator, i.e., terms proportional to S% and Sy. The 
lowest order transverse terms must be proportional to 
Si_ + S*i.LJ If these were the only terms in the Hamil- 
tonian responsible for tunneling, then only resonances 
with "tunneling length" m! — m (see Fig. |l|), which is 
a multiple of four, would have been observed. This is 
not the case for Mni2. To explain the presence of res- 
onances with all k, one should invoke a transverse field 
contribution to the Hamiltonian, —H X S Z . Dipolar fields 
from magnetic moments of Mni2 molecules and hyper- 
fine fields fr^TLjMjt^pcLai .have been suggested as natural 
candidates .oEjuj'OCjOO Experiments, however, indi- 
cate that transverse fields needed to explain the data are 
stronger than the ones produced by atomic nuclei and 
magnetic dipoles. The nature of the effects responsible 
for tunneling in Mni2 remains an open question. 

Next controversy is the origin of the \ft law that ap- 
proximately fit^tbje initial stage of the relaxation in some 
experimentsBUS Suggested explanations include col- 
lective efEects due to dipolar interaction .between Mni2 
moleculesEi and fluctuating random noiseO The first re- 
quires some specialpinitiai conditions which are not sat- 
isfied in experimentOO'Cll while the second requires cer- 
tain assumptions about the spectrum of fluctuations. 

Another open question is the nature of the "minor 
species" of Mn.12 that relaxes faster than the "major 
species", and is present in all samples studied to date. 
In the subkelvin temperature range, when the relaxation 



slows down, the minor species is the main source of the re- 
laxation seen in experiment. Wernsdorfer et al attribute 
that species to the defective sites of the crystal lattice 
(see Footnote 9 in Ref. ^). On general grounds, it is 
obvious that defects inside the Mni2 clusters cannot ac- 
count for the minor species. Such defects would strongly 
modify the spin Hamiltonian, which would result in a 
significant change of the resonant fields and orders-of- 
magnitude change of tunneling rates, in contrast with 
experimental observations. 

On the contrary, defects of the tetragonal crystal lat- 
tice such as dislocations, which do not change the struc- 
ture of the cluster, should produce symmetry-violating 
terms in the Hamiltonian that are responsible for tunnel- 
ing. Some of these terms are quadratic in S x and S y . For 
even k [see Eq. (^)] they produce tunneling in a lower or- 
der of the perturbation theory than the transverse-field. 
In addition, dislocations give rise to local rotations of the 
easy axis, which for H z results in a transverse field 
that unfreezes resonances with odd k. 

Because of the long-range nature of deformations 
caused by dislocations, the number of affected crystal 
sites should be relatively large even for a moderate con- 
centration of dislocations. Strong deformations will exist 
only at a small number of sites inside dislocation cores. 
Most of the Mni2 molecules will develop weak deforma- 
tions, with the tunneling rate depending on the location 
of the molecule. Consequently, the relaxation process 
evolves from the relaxation of the minor species, close 
to the dislocation cores, to the relaxation of the major 
species far from the dislocation cores. 

In this paper we show that edge and screw dislocations 
(see Fig. |^) should be the main source of spin tunneling in 
Mni2 crystals. Broad distribution of deformations causes 
broad distribution of thermal activation and tunneling 
rates. To study the magnetic relaxation in Mni2, we de- 
velop the theory of incoherent Landau-Zener tunneling 
from excited energy levels. We compute the relaxation 
law in a crystal with dislocations and show that it obeys 
scaling that can be tested in experiment. We also demon- 
strate that in the kelvin temperature range the relaxation 
is well described by the stretched-exponential law. 

The paper is organized as follows. Spin-lattice cou- 
plings due to different types of dislocations are studied 
in Sec. [ij]. Sec. [II contains perturbation formulas for the 
tunneling level splittings for even and odd resonances. 
The analytical formula for the distribution of transverse 
anisotropics due to a random array of dislocations is de- 
rived in Sec. IV. The distribution of tunnel splittings in 
Mni2 crystals is studied in Sec. Sec. [Vl| is devoted 
to the incoherent Landau-Zener transitions at finite tem- 
peratures. In Sec. VII the results are applied to Mni2 



cryst als with the distribution of tunel splittings. In Sec. 
VIII we compute the magnetic relaxation in the ther- 
mal activation regime. The implications of the results 



obtained in this paper are summarized in Sec. IX 
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FIG. 2. Edge and screw dislocations. 



SPIN-LATTICE COUPLINGS DUE TO 
DISLOCATIONS 



A. Spin Hamiltonian 



We study the Hamiltonian 

H = -DSl - H Z S Z 



TLn 



(2) 



where S z is the z-component of the spin operator, S = 10, 
D = 0.65 K, H z is the magnetic field applied along the 
z-axis (c-axis of the crystal), and 7i mc is the magnetoe- 
lastic coupling. The Hamiltonian of Mni2 also contains 
crystal fields of fourth order on the spin operator, mag- 
netic dipole interactions, and hyperfine interactions. We 
neglect them in order to emphasize the effect of dislo- 
caticoiK|_X]jie magnetoelastic coupling in Mni2 is of the 



H mo = giD(e 

xx ^yy )( s x - s v ) + 92De xy {S x , Sy} 



D(a x {S x , S z } + a y {S y , S z }), 



where {A, B} is anticommutator, 



33 e. 



X V — 93 £ yz + 9A<JJyz, 



and 



1 / du a 
= 2 \dxp~ 



dx a 



Wa0 = 



1 / du a 
= 2 



du fj 

8x a 



(3) 



(4) 



(5) 



are linear deformation tensors; oj a p being a pure rotation 
and u being the displacement. Note that, in principle, 
the Hamiltonian may contain terms of higher order on 
deformations. The exact strain tensor also has quadratic 
contribution on displacement. These terms may become 
important at and near dislocation cores, that is, for a 
small number of Mni2 molecules which do not contribute 
much to the relaxation process. 

The effect of dislocations is determined by values of 
constants g\, <?2j<?3j an d <?4- The.Kahie of is fixed by the 
symmetry of the crystal, g± = IoIIjO There are general 
arguments that all ^-constants must be of the same order, 
see Ref. ^0] and references therein. For numerical work 
we will use g\ — gi = .93 = g± = 1. 



Before calculating tunneling rates of Mni2 molecules 
due to dislocations, it is convenient to make rotations 
of the spin axes in (x,y), (x,z), and (y, z) planes in 
order to transform the Hamiltonian of Eq. (|2|) into a 
diagonal quadratic form on spin operators. The terms 
{S X ,S Z } and {Sy,S y } in Eq. (|^) are eliminated by ro- 
tations by small angles, a x /2 and ot y /2, in the {x, z) 
and (y, z) planes, respectively. This results in a slight 
renormalization of the uniaxial anisotropy constant D 
and of the transverse-anisotropy terms. Both effects are 
quadratic on deformations and are neglected in our ap- 
proach. The rotations also change the form of the Zee- 
man term —H Z S Z , resulting in the transverse field 



H± = 



1 



a 2 y H z 



(6) 



in the rotated coordinate system, which is linear on de- 
formations and proportional to the magnetic field H z ap- 
plied along the c-axis of the crystal. This transverse field 
unfreezes tunneling transitions which change the spin 
projection m by an odd number. The terms {S x , S y } in 
Eq. (Q) are eliminated by rotations in the (x, y) plane by 
an angle which is not small. They produce the transverse 
anisotropy E(S^,—Sy,) in the first order on deformations. 

In the four subsections below, we will consider differ- 
ent types of linear dislocations. We will use formulas 
of isotropic elastic theory for displacements,^!! which are 
accurate enough for our conclusions and, at the same 
time, are less cumbersome than the exact expressions for 
the tetragonal crystal symmetry. We will see that for 
all types of linear dislocations the spatial dependence of 
the transverse anisotropy and the transverse field can be 
written as 



E = 2D 



H± = H : 



(7) 



where r is the distance from the dislocation axis, mea- 
sured in the lattice units, whereas g(cp) and 3j?(<^) are 
functions of the angle, which are of order one if gi ~ 
92 ~ 93 ~ 94 = 1< One can immediately see from Eq. 
(0) that the effect of dislocations on tunneling must be 
strong. Indeed, for r ~ 1 one has E ~ D, whereas the 
spatial decay of E is slow, so that each dislocation affects 
a large number of Mni2 molecules in the crystal. 



B. Screw dislocation along the anisotropy axis 

For screw dislocations shown in Fig. [2] the only nonzero 
component of the displacement u is m z . For the disloca- 
tion axis at x, y — it is given by 



, . c y 
Uz{x 7 y) = — arctan-, 
Zir x 



(8) 



where c is the lattice spacing along the c-axis. The only 
nonzero components of deformations are 



y 



4ir x 2 + y 2 

c x 

vz 47r x 2 + y 2 



(9) 



■1 



Thus, according to Eq. (|3J), screw dislocations along the 
c-axis do not produce transverse anisotropy of the first 
order in deformations. They, nevertheless, create a trans- 
verse field in the rotated coordinate system, which is 
given by 



C |g3 — g4 1 

8tt yj x 2 + y 2 



H,. 



(10) 



[cf. Eq. @]. This transverse field gives rise to spin 
tunneling of order H™ ~ m between the resonant pair 
of levels e m and e TO < [see Eq. (2.6) of Rcf. Since 
m'—m is a usually a large number, this tunneling is much 
weaker than tunneling due to the transverse anisotropy 
E from other types of dislocations, which appears in the 
£>(m -m)/2 orc j er f or even resonances (see Sec. fill for more 



detail). Since all types of dislocations must be simultane- 
ously present in a crystal, we conclude that c-axis screw 
dislocations are irrelevant for spin tunneling. 



C. Screw dislocation perpendicular to the anisotropy 

axis 

For a screw dislocation along the Y-axis the only 
nonzero component of the displacement is 



u v (x.z) = — arctan — 



(11) 



where b is the appropriate lattice spacing. Thus the 
nonzero components of the deformation tensors are 



x 



" " 47T X 2 + Z 2 

b Z 



--xy 



An x 2 + z 2 



(12) 



To bring the spin Hamiltonian to the canonical form, one 
should first rotate the spin axes in the (y, z) plane by an 
angle ot y /2, which results in a transverse field 

H± = \a y H z = (.93 + 9i)^-^-^H z (13) 
2 87T x z + 

directed along the new Y-axis. Then the rotation by <f> = 
7r/4 in the x,y plane transforms {S x ,S y } into —(S 2 , — 
S 2 ,), which results in the Hamiltonian 

TL = —DS 2 Z , — H Z S Z > + E(S X , — S y i) — H x rS x > — H y /S y i, 

(14) 



where 



and 



E = - 92 De xy = a ,D-^— i 



H x > = H y , = (g 3 + 54) 



8V2n x 2 + y 2 



(15) 



(16) 



Transverse anisotropy: Screw Dislocation j/Y 




5 10 15 

x, lattice units 



FIG. 3. Contour plot of transverse anisotropy E created 
by a screw dislocation along the Y-axis, see Eq. ©. (Grey 
scales are arbitrary.) 



(the signs of H x i and H y i are irrelevant). This type 
of dislocations generates the spatially-dependent trans- 
verse anisotropy which is the main source of tunneling. 
It also generates a transverse field that unfreezes odd 
resonances. The hard and medium axes, x' and y' , in- 
terchange when one crosses the plane 2 = 0, whereas the 
transverse field always has equal components along the 
hard and medium axes. 



D. Edge dislocation along the anisotropy axis 

For this dislocation shown in Fig. || there is no dis- 
placement along the c axis, i.e., u z = 0, whereas other 
displacement components are given by 



b 

2^ 



arctan ■ 



xy 



2(1 - a) x 2 +y 2 



(17) 



and 



b 

"2^ 



1 -2a 

4(1 -t) 



ln(x 2 + y 2 ) + 



1 



2(1 - cr) x 2 +y 2 



(18) 



where 0<er<l/2is the Poisson elastic coefficient (we 
will use a — 0.25 in the numerical work). For the defor- 
mations one obtains 



x 2 y 



and 



-xy 



7r(l - a) (x 2 + y 2 ) 2 

b x(x 2 — y 2 ) 
'4tt(1-ct) (x 2 +y 2 ) 2 ' 



(19) 



(20) 



It is clear that this dislocation does not tilt the anisotropy 
axis and thus does not generate a transverse field. Af- 
ter the rotation in the (x, y) plane the spin Hamiltonian 
becomes 



,1 



Transverse anisotropy: Edge Dislocation \\Y 
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FIG. 4. Contour plots of transverse anisotropy E created 
by: a - an edge dislocation along the F-axis, see Eq. (p8[); b - 
an edge dislocation along the Z-axis, see Eq. 



T~i — —DS Z — H Z S Z + E(S x r — Sy') 



with 



£ = DJgf(e xx - e yv ^ ^ ~ 2 - 2 



92 £ xy 



(21) 



(22) 



E. Edge dislocation perpendicular to the anisotropy 

axis 

For such dislocations the nonzero displacements are u x 
and u z . They are given by Eqs. ( p~7|) and (|lj) with y => z. 
For deformations one obtains = 0, 



/. J3 - 2q)z 2 + (1 - 2a)z 2 



47r(l - a 



{x 2 + z 2 ) 2 



and 



/ . ,du x du z 
(33 + 94) -Q- + (93 - .94) -7^7 



(24) 



"57 ~~ 4tt(1-ct) : 



and 



<9u 2 



(3-2a)a; 2 + (l-2o-)z 2 
(x 2 + z 2 ) 2 

(1 - 2a)x 2 + (3~2a)z 2 
~dx~ ~ 4Tr(l-a) U ' (a; 2 + z 2 ) 2 ' 

The resulting spin Hamiltonian has the form 

H = —DS Z , — H z S z i + E(S X , — Sy,) — H x iS x 

where 



E = gi D(e a 



-a x H z . 



(25) 



(26) 



(27) 



(28) 



Above the XY plane (at z > 0) the transverse field is 
directed along the hard axis, while at z < the transverse 
field is along the medium axis. 



III. TUNNEL SPLITTINGS 

For concentrations of dislocations upto c = 10~ 2 (in 
lattice units) most of Mni2 molecules are still far from 
the dislocation cores. Thus the transverse anisotropy and 
the transverse field generated by dislocations [see, e.g., 
Eq. (p7|)] are typically small in comparison with the uni- 
axial term in the Hamiltonian of the ideal crystal. In this 
situation one can obtain the tunneling splittings A mm / 
of resonant pairs of levels in the lowest-order of the per- 
turbation theory.. Epr, . th e model with the transverse field 



the splittings are 



2D 



[(m' — m — l)!] 2 

1 {S + m>)\{S - m)\ 
(S - m')\(S + m)\ 



H x 
2D 



(29) 



For the transverse anisotropy model, the splittings have 
been calculated in Rcf. |53| for H z — and generalized in 
Ref. ^ for any field bias. It is convenient to rewrite Eq. 
(4) of Ref. |(] in the form 

2D 



[0 



2)!!] 2 



' (S + m'y.jS-my f\E\\ (m '~'" ! ' 2 
(S - m')l(S + m)\ \2D, 



(30) 



In the transverse-anisotropy model, A mm i is nonzero for 
even m! — m = —2m — k, which for even S requires an 
even value of the tunneling resonance number k. For 
odd values of m! — m, tunneling resonances are quenched, 

In Mni2 samples with dislocations, both the transverse 
field and the transverse anisotropy are present, see Eq. 
(0). Since at resonances the longitudinal field satisfies 
H z = kD, the transverse field is well defined for a given 
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FIG. 5. Exact and perturbative results for the trans- 
verse-field dependence of the splitting of the metastable 
ground-state at the k = 7 tunneling resonance. 



type and configuration of dislocations, with E and H± 
being of the same order of magnitude. To obtain tunnel 
splittings perturbatively in such a model, one has to sum 
up the products of matrix elements of perturbations and 
corresponding energy denominators, each product—rep- 
resenting a chain connecting the states to and m'.LJ In 
general, this algorithm is very cumbersome and has no 
advantages with respect to direct numerical diagonaliza- 
tion of the spin Hamiltonian. 

Fortunately, as will become apparent below, for H± ~ 
E <C D the effect of the transverse field on tunneling 
is much weaker than that of the transverse anisotropy. 
Thus, for the even resonances one can neglect the field 
contribution to the splittings, whereas for the odd reso- 
nances the field should be taken into account in the first 
order only. In the latter case, the main source of the split- 
ting is the transverse anisotropy taken in a high order of 
a perturbation theory needed to build a chain of matrix 
elements that changes the spin projection by ml — m— 1. 
For odd resonances the transverse field makes the miss- 
ing single perturbation step along that chain, needed to 
change the spin projection by one. The corresponding 
matrix element (mi\H x S x \mi + 1) can be inserted at any 
place in the chain, mi = m, m + 1, . . . , m! — 1. Quite 
fortunately, the corresponding sum can be calculated ex- 
actly: 



1 - 



p=0 



r(r-p+l/2)r(p+l/2) 

r(r-p + i)r(p + i) 



= i 



(31) 



and the resulting splittings for odd resonances acquires 
the elegant form: 



H , 



[(m' 



2)!!] 2 



>'V(5- 



r'W 



\E\_\ 
2D J 



(m 1 — rn — 1) /2 



(32) 



Note that this perturbative result does not depend on 
the sign of E, that is, the splitting is the same for the 



transverse field directed along the medium or hard axes. 
The difference between these two cases appears in higher 
orders in the transverse field, see Fig. ||. For the trans- 
verse field along the hard axis, there is a destructive in- 
terference of the contributions of the transverse field and 
the transverse anisotropjO that cause splittings to van- 
ish atE3'E3 



H x = k x ^2E(D + E). 



(33) 



For the integer spin S the number k x satisfying the con- 
dition \k x \ + \k\ < 2S + 1 is odd, where k is the tunnel- 
ing resonance number of Eq. ([!]) . For even resonances k 
splittings vanish at k x — 0, ±2, . . ., whereas for odd reso- 
nances the splittings vanish at k x = ±1, ±3, ... It is clear 
from Eq. (33) and from Fig. ||that the linear dependence 
of the splitting on the transverse field, as given by Eq. 
(|32|), is valid for H± <C V2ED, which is satisfied in our 
case Hj_ — E < D. 

The condition on the transverse field obtained above 
also follows from simple perturbative arguments. For 
the matrix element of the perturbation V mi:mi +2 = 
(mi\V\mi + 2) one has V mi j ni +2 ~ E if the perturba- 
tion is due to the transverse anisotropy, and V r mii „ ll +2 ~ 
H^_/D if it is due to the transverse field. Thus the 
effects of both types of perturbation are comparable if 
H x - VED. 

Let us compare now the strengths of even and odd 
tunneling resonances. For an odd resonance the ratio of 
A mm ' to the geometrical mean value of the splittings of 
adjacent even resonances has the form 



A, 



F(S, TO, TO,') 



H,, 



V2ED' 



\f ^Tn,m' ,rn' — l 

where 

(to' — to — 1)!!(to' — to — 3)!! 



(34) 



F(S, to, to') 



[(to' -to- 2)!!] 



(S + m')(S - m! + 1) 



{S-m')(S + m' + 1) 



1/4 



(35) 



The function F(S, to, to') is of order unity. It approaches 
7r/2 in the quasiclassical case of large S, to, and to'. Thus 
the odd resonances are weaker than the even ones by the 
factor of order H x j\]DE. This effect is not dramatic 
since odd and even resonances only differ by a prefactor 
in front of high power of the small ratio E/D, see Eqs. 
© and @. 

It should be stressed that for numerical calculation of 
tunnel splittings the exact resonance conditionE3E3 



H z = k\/D 2 - E 2 , jfc = 0, ±1, ±2, . 



(36) 



should be used instead of Eq. |T]). Although for E <C D 
the corrections to the resonance fields H z due to the 
transverse anisotropy are small, the detuning of levels 
arising from the use of the approximate resonance con- 
dition of Eq. (|) is much larger than the level splitting. 
If the exact Eq. (pq) was not known (which may well be 



the case for more complicated spin models), the numer- 
ical finding of level splittings would require sweeping of 
H z which is a time-costly procedure. 

One of the results of this section is that S > 1 the 
slowly decaying perturbations of the spin Hamiltonian 
due to dislocations, Eq. (|?]), generate tunnel splittings 
which change by many orders of magnitude on the dis- 
tance from the dislocation. Consequently, the behavior 
of Mni2 crystals at low temperatures must be similar to 
that of disordered systems with widely distributed pa- 
rameters. The natural scale to discuss physical prop- 
erties of such systems is logarithmic, which makes the 
results practically insensitive to any prefactors of order 
unity, such as the functions g(cp) in Eq. (Q) encapsulating 
differences between different types of linear dislocations. 
For this reason, we will use only one type of dislocations 
for illustrations, the edge dislocations along the y-axis. 
Before studying spin-tunneling rates due to dislocations 
in real Mni2 crystals, one should work out one funda- 
mental characteristic of these crystals: the spatial distri- 
bution of transverse anisotropies. This will be done in 
the next section. 



IV. DISTRIBUTION OF TRANSVERSE 
ANISOTROPIES DUE TO DISLOCATIONS 

In a crystal with dislocations, the deformation tensor 
at any given point is a sum of contributions due to many 
different dislocations. The superposition principle for de- 
formations follows from-the linearity of the equations of 
the theory of elasticity^ and it holds everywhere out- 
side dislocation cores, i.e., for the distances from the 
dislocation axes r > 1. Statistical properties of defor- 
mations and thus of the spin tunneling rates in a Mni2 
crystal depend on the spatial distribution of dislocations 
which is poorly known. In this section we present analyt- 
ical solution of the problem of distribution of transverse 
anisotropies, assuming that dislocations are distributed 
at random. Note that alternatively, one can consider a 
regular array of dislocations with alternating directions 
of the Burgers vector, to achieve a balance of contractions 
and dilations throughout the crystal. Such a dislocation 
array should be randomized to some extent to make it 
more realistic, as was done in Ref. [Tl[ 

For the random-dislocation model, one should, in prin- 
ciple, consider different types of dislocations and compute 
an average over all crystal sites. For simplicity, we will 
take into account only one type of randomly distributed 
dislocations with collinear axes: the edge dislocations 
along the F-axis, see Sec. [n]E, and consider one repre- 
sentative site in the middle of the crystal, r = 0. The 
distribution function for the transverse anisotropy can 
be defined as 
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2D' 



(37) 



where N ^> 1 is the number of dislocations in the crystal 
and the averaging is carried out over their positions in 



the plane perpendicular to the dislocation axis within a 
circular region of radius R. One can define 
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where c is the concentration of dislocations and R c is the 
characteristic distance between dislocations. 

Let us at first analyze the large- \E\ asymptotes of 
due to the regions with large deformations of both signs 
close to one of dislocations. In that case one can neglect 
the influence of all other dislocations and consider the 
one-dislocation model 
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Integration yields 
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(40) 



where E c is the characteristic transverse anisotropy at 
the distance R c . This formula becomes invalid for E < 
E c , where the lines of constant E in Eq. (|3^) cross the 
boundary of the region under consideration, r = R c . In 
fact, for E < E c the very Eq. ( |39"| ) becomes invalid and 
one has to take into account other dislocations. Eq. ( |40| ) 
suggests that one should introduce the distribution func- 
tion for a reduced transverse anisotropy a 



fa = E c fg, a = E/E c , 
which has the asymptote 

fa = 1 joi a ^ 1- 

\a\ 6 



(41) 



(42) 



In the general case, with the help of the identity 
2ttS(x) = duie luJX , the averaging over the coordinates 
of different dislocations in Eq. fl37|) can be factorized, 



J L = I ^e^f(u) N , 



where 
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As we shall see, in Eqs. (g|) and Q), to ~ R c < R for 
N 3> 1, thus the argument of the exponential in Eq. (|44|) 
is small and f(u>) is close to unity. Then the exponential 
can be expanded and integrated, with a log accuracy, in 
the interval \u>\ < r < R. Given that (g(<p)) = 0, the 
result has the form 
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where Co is a constant of order unity. Now with the use 
of Eqs. (|3^) and (E3) one can write 
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At this point one may forget about the initial assumption 
on the circular form of the spatial region. The shape of 
the crystal only affects the value of the constant Co under 
the logarithm. Eq. (Q) confirms the assumption lu ~ 
1/E C ~ R c made above. Now we are prepared to write 
down the final result which is convenient to formulate in 
terms of the function f a defined by Eq. (fid]) 
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Here the cutoff A satisfies 1 C A C vAf; one cannot 
integrate up to oo since the form if the integrand is only 
valid for u <C vN. Clearly, for large enough crystals with 
N 1 the result does not depend on A. We remind that 
for the edge dislocations along the F-axis, the distribu- 
tion of transverse anisotropics is an even function. This 
distribution is shown for E > in Fig. |(J. 

Integrating Eq. ( f47|) by parts three times, one can re- 
cover the asymptote of f a at \a\ S> 1 which is given by 
Eq. (^). This power-law asymptote is a consequence of 
the logarithmic singularity of the integrand in Eq. ( [47] ) 
at u — > and it leads to the divergence of the second 
moment of f a . On the other hand, for large N the distri- 
bution function may be well approximated by Gaussian 
for not too large a. Indeed, for large N the logarithm in 
Eq. (47) is weakly dependent on u and can be replaced 
by a constant. The best value of this constant corre- 
sponds to u for which the argument of the exponential 
equals one. This requires solving a transcendental equa- 
tion that can be done in a perturbative way. With a good 
accuracy one can use 
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which results in the aproximation 



fa = 
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: exp 
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(48) 
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which is also shown in Fig. ^|. 

As the number N of dislocations in the crystal in- 
reases, the function f a of Eq. (|47]) becomes closer and 
closer to the Gaussian, whereas the power-law asymptote 
given by Eq. (^|) becomes shifted to the region of very 
large a where it is hardly visible. This effect is due to 
the accumulation of small contributions from dislocations 
situated at large distances from the observation point 
(of order of the linear dimension of the crystal). Such 
small contributions from distant dislocations, which lead 
to the Gaussian distribution f a , win over contributions 
from close dislocations responsible for Eq. (ff2|). One 
can estimate the characteristic value a\ of the transverse 
anisotropy at which the distribution function changes its 
behavior by equating the Gaussian approximation for f a 
given by Eq. (fl|) to its asymptote of Eq. (H|). This 
yields 
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FIG. 6. Distribution of the transverse anisotropy due to a 
random array of edge dislocations along the Y-axis 



ai = J 4Zln(4L/\/7r)- 



(50) 



Transverse anisotropics E [see Eqs. (|37]) and ([II])] cor- 
responding to a > a± are due to a single dislocation in 
the vicinity of the observation point, whereas a < oi\ 
are collective contributions of many distant dislocations. 
The characteristic distance r \ from t he dislocation core 
is defined by E\ = a\E c = y/ (g(tp) 2 ) jr\ and reads 
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(51) 



One can see that in a macroscopic crystal the contribu- 
tion of a particular dislocation dominates within a dis- 
tance from its core, n, that is small compared to the 
average distance between the dislocations. The fraction 
of cites in the lattice affected mostly by one close dislo- 
cation is 
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(52) 



This fraction is small if the number of dislocations N in 
the crystal is large. Molecules belonging to this group 
may be interpreted as the minor species. One should 
note that the value of n\ above is, in fact, only the upper 
bound on n\ . The asymptote f a = I /a 3 is applicable for 
a 3> oi\ (say, for a ^ 3ai) rather than for a ^ ai. If 
one replaces ot\ by 3ai in Eq. (52), the value of n\ will 
decrease by one order of magnitude. 

The Gaussian approximation for the function j ^ with 
the help of Eq. (|4l|) can be written in the form 
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where 
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That is, the accumulation of contributions from distant 
dislocations leads to the effective logarithmic renormal- 
ization of the concentration of dislocations c with L de- 
fined by Eq. (ES) . For edge dislocations along the Y-axis, 
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FIG. 7. Distribution of the logarithm of transverse 
anisotropy due to random array of edge dislocations along 
the Y-axis 



the quantity y/(g((fi) 2 ] 
( p3| ) is given by 



according to Eqs. (0), (H|), and 
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where a is the Poisson elastic coefficient. For for g\ = 1 
and a = 0.25 one has ^f{g{tp) 2 ) « 0.084. 

The experimentally studied Mni2 crystals are rather 
large, about 0.5 x 0.5 mm 2 , which corresponds to the 
cross-section of about 10 11 lattice cells. Even for the con- 
centration of dislocations as small as c = 10 -4 per cell, 
the number of dislocation in the crystal is about N w 10 7 . 
For co = 1 this gives L — 9.1, i.e., the effective concentra- 
tion of dislocations increases by an order of magnitude, 
c = 0.91 x 10~ 3 . The corresponding value of E £ that fol- 
lows from Eqs. (fj) and © is E- c = 0.449 x 10~ 2 . For 
c = 10~ 3 one obtains L = 10.3, thus E e = 1.51 x 10~ 2 . 
If one takes L m 10, then Eq. @ yields n x w 1/125. 
That is, less than 1% of all Mni2 molecules in the crys- 
tal belong to the fast relaxing group in the vicinity of 
a single dislocation. In fact, the value of m should be 
much smaller, see comment after Eq. (|52"|). The renor- 
malization of the concentration of dislocations and the 
Gaussian distribution of transverse anisotropics for large 
crystals are clearly seen in Fig. ^: The distribution be- 
comes broader in the ct-scale due to the increase L with 
N, Eq. ©. 

Having determined the distribution of transverse 
anisotropies created by dislocations in a Mni2 crystal, 
we can calculate the distribution of tunnel splittings for 
different level pairs and different resonances. This will 
be done in the next section, where we will use the distri- 
bution of the logarithm of the transverse anisotropy 
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f L (x)=2e x fE(e x ), x = \nE. 



If the Gaussian approximation for /g, Eq. (|53[), is 
adopted, then fa has the form 
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This function is plotted in Fig. for c = 10 3 (Ec 
1.51 x 10~ 2 ) and c = 10" 4 [E~ c « 0.449 x 10~ 2 ). 



V. DISTRIBUTION OF TUNNEL SPLITTINGS 

Distribution of transverse anisotropies due to disloca- 
tions, which was obtained in the previous section, deter- 
mines the distribution of the tunnel splittings which for 
the even-fc resonances are given by Eq. ( |30| ) . It is conve- 
nient to rewrite Eq. ( |30| ) in terms of the resonance num- 
ber k and the level number n — 0, 1, ... in the metastable 
well of negative spin projections 

m = n — S, m = S — n — k. (58) 

Thus Eq. (|3C]) becomes 

A nk =g nk E^ k (k even), (59) 

where E = E/(2D), £ nk = S - n - fc/2, and 
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With the help of Eq. ( J5£| ) one can write down the distri- 
bution of the square of tunnel splittings for any resonance 
in terms of the previously defined 
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(61) 



Because is a large number, especially for resonances 
with small n and k, tunnel splittings in a Mni2 crystal 
spead over many orders of magnitude, from rather large 
values near dislocation cores to very small values far from 
dislocations. It is thus more convenient to consider the 
distribution of the decimal logarithm of tunnel splittings, 
y = log 10 A. With the help of Eq. (M) it can be written 
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(62) 



Here Jl (x) is the distribution function of x = In E 
defined by Eq. (p6|). Distribution of the ground-state 
(n = 0) tunnel splittings for different resonance numbers 
k is shown in Fig. H for a Mni2 crystal of size 0.5 x 0.5 



mm with concentration of dislocations c = 10 



One 



can see that the distribution is shifted to the left and 
becomes broader for smaller k. 

For odd tunneling resonances, the splitting, according 
to Eq. (32), depends on both the transverse anisotropy 
and the transverse field which are generated by disloca- 
tions. Since the functions g(<p) and gnif) hi Eq- (0) are 
different, one cannot, strictly speaking, express the distri- 
bution of splittings via that of the transverse anisotropy 
alone. Fortunately, transverse field is only a prefactor 
in front of a small term in Eq. (p3). This allows one to 
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express H± via E with the help of Eq. (0) and use the 
resonance condition H z = kD to obtain 
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with the ratio of angular functions approximated by a 
constant c oc jd ~ 1. In this approximation, the level split- 
tings at odd k are functions of the transverse anisotropy 
alone, so that one can rewrite Eq. as 



A nk = g' nk E^ 



where 
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[cf. Eqs. (|59| ) and (p0|)1. For the distribution of the deci- 
mal logarithm of splittings at odd resonances one obtains 
the formula analogous to Eq. (^) where the constant c 
enters under the logarithm. Note that the power £' nk for 
an odd resonance is the same as the power for the 
preceding even resonance, thus the splitting distribution 
functions for these two resonances differ, according to Eq. 
(|62|), only by a shift. This is clearly seen on Fig. || where 
we have plotted the odd-resonance curves for c Q dd = 1. 

If the distribution of the level splittings for different k 
is simulated, as was done in Ref. |d], or is measured, the 
theory developed above could be tested by re-plotting the 
data in terms of the scaling variable^ 
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Since x = In E, it does not depend on the numbers n and 
k, thus all the curves should scale. For odd resonances, 
one should use £' nk and g' nk of Eq. (|6|) , where c can be 
considered as an adjustable parameter. 



VI. INCOHERENT LANDAU-ZENER 
PROCESSES 

A very convenient method to experimentally study 
spin tunneling consists of sweeping the longitudinal field, 
thus making a pair of levels on different sides ofrtjie bar- 
rier, e n and £„', to go through a resonance. ErEZI For a 
purely quantum-mechanical system, the probability of 
the system to stay in the same well (i.e., to go „ fr fm one 
exact energy branch to the other) is given by 
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exp 
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(67) 



where A„ is the level splitting and v n = \d(e n > — e n )/dt\ 
is the speed of the level detuning unperturbed by the tua-, 
neling interaction, i.e., the energy sweep rate. LandauL3 
has obtained this formula in the special case of a large 
argument of the exponential (low sweep rates), whereas 
ZenerE3 obtained it for all sweep rates, however with a 
wrong additional factor 2ir in the exponential. Recent 
Ref. Pq reproduces Zener's results without the wrong fac- 
tor 2tt [Eq. (11) of Ref. ||is equivalent to Eq. @]. We 
will call this process the coherent Landau-Zener process. 

For systems interacting with the environment, the sit- 
uation changes if the linewidth T n of the level is greater 
than the level splitting A„. This, in fact, happens in 
most cases, since A„ for spin systems is a high power of 
a small perturbation, while T n is due to transitions be- 
tween the levels in the same well which appear in much 
lower order of the perturbation theory on spin-phonon 
(or other) interactions. The only exception is the tunnel- 
ing resonance between the ground-state levels in the two 
wells at low temperatures because the phonon processes 
contributing to their linewidths die out exponentially at 
T — > 0. If at least one of the resonant levels is not the 
lowest level in the well, it has a considerable linewidth 
down to T = due to the transitions onto lower levels in 
the same well, accompanied by the emission of phonons. 
In this case T„ 3> A„ and tunneling becomes mediated 
by the environment, i.e., incoherent. 

Incoherent spin tunneling should be described in terms 
of the density matrix emiation rather than in terms of the 
Schrodinger equation.Ej Since at low temperatures the 
escape process (via tunneling or via thermal activation) is 
much slower than the equilibration within the wells, one 
can eliminate the non-diagonal elements of the density- 
matrix equation and obtain the system of equations for 
the level populations that describes tunneling, Eq. (4.12) 
of Ref. [l^. In sweeping-field experiments on Mni2 in 
the kelvin or subkelvin range, the backflow of particles 
from the stable to the metastable well is exponentially 
small and can be safely neglected for all resonances with 
k 7^ 0, since the corresponding activation energy is about 
2SD ~ 13K. Thus there is only one-way escape from the 
metastable well described by the equatior£3 
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r„/ , the numbering of levels begins 



from the metastable ground state, n 



S, and n top 



i n 



corresponds to the top of the barrier. The tunnel split- 
ting A n k depends on both the level number n in the 
metastable well and the number of the resonance k of Eq. 
(HI). They are given by Eqs. (|Fj) and @ with m = n-S 



and 



= S- 



k. In the kelvin and subkelvin temper- 



ature range almost all Mn 12 molecules in the metastable 
well are in the ground state, whereas the populations 
of the excited states entering Eq. ( |68|) are exponentially 
small and are given by the equilibrium formulas 



N n = N e 
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where the level energies are taken at the resonant values 
of the field, H z = kD (k > 0) 



e n = Dn(2S - n - k) - DS(S - k). 



(70) 

Now Eq. d6§| ) can be integrated for the field sweeping 
across the resonance. For the constant energy-sweep rate 

£«' - = v nk t + const, (71) 

where v n k = (2s — 2n — k)dH z /dt, one obtains 
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where _/y ( ( fe > boforc ) arLC j jy | ( fc > after ) are numbers of Mni2 
molecules in the metastable well before and after crossing 
the fc-resonance. 

Eq. (72) describes a superposition of incoherent 



Landau-Zener tunneling processes which come both di- 
rectly from the metastable ground state and via excited 
states. At T = the result coincides with that of Eq. 
(|67|), although the physics is different. Note that for 
the incoherent Landau-Zener tunneling between the two 
lowest-energy states one has to take into account the flow 
of particles in both directions, which results in the disap=. 
pearance of the factor 2 in the denominator of Eq. (f72|) .c2l 
In this situation, however, one should carefully check the 
applicability condition of the method, F„ 3> A n , since 
r„ may be very small. Incidentally, for Mni2 acetate 
the tunneling between the ground states is so weak that 
the k = resonance cannot be observed experimentally, 
unless a strong transverse field is applied. 

The field sweeping technique is a very convenient ex- 
perimental tool since Eq. (72) does not depend on the 
damping parameters which are difficult to determine with 
sufficient accuracy. The combination A^ fe e~ Ae " fc / T is fa- 
miliar from Refs. 13, ^FJ, ^l] where it was used to inves- 
tigate the transition on temperature between thermally 
activated and tunneling regimes of the escape from the 
metastable well. To be more precise, let us rewrite Eq. 
(|72|) with the help of 
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In fact, this form of f nk is valid when the tunnel splitting 
is small, A nk -C Se„k ■ This is true for energy levels below 
the top of the barrier. Near the top of the barrier, where 
A„fc ~ 5e n k, a more accurate formula isc3 
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It ensures that /„j. ~ Ssn^e" nk near the top of the 
barrier for any approximation used for A n fc. 

Since f(n) is a product of the two functions, one of 
which rapidly increases and the other rapidly decreases 
as n goes up, there is a narrow group of levels around 
n = n*(T) that maximizes f n k(T) and makes the domi- 
nant contribution to the Landau-Zener transition. Above 
the quantum-classical transition temperature To, one has 
n*(T) = 7i t op, which corresponds to the activated regime 
with the spin escaping over the top of the barrier. As T 
is lowered below To the tunneling level goes down, which 
corresponds to the thermally assisted tunneling or, at 
even lower temperatures, to the ground-state tunneling. 
If n*(T) is continuous, one can speak about the second- 
order quantum-classical transition. If n*(T) jumps over 
several levels at some temperature, the transition is first 
order. |l3|, ^FJ, |2l] Recently an experimental evidence of 
the first-order transition in Mni2 has been obtained in 
Refs. [n], [12]. (The biaxial model with the magnetic field 
along the hard axis shows_two first-order transitions in 
some range of parameters.Ej'EJ) 

Experimental study of the temperature dependence of 
the tunneling level n*(T) in Mni2 is possible due to the 
small quartic-anisotropy term — AS Z in the spin Hamil- 
tonian, which makes the resonance fields H z dependent 
of the quantum number n. Tunneling resonance between 
higher pairs of levels takes place at slightly lower value 
of H z than the resonance between lower pairs of levels, 
that is, the fc-resonance has a fine structure. Measuring 
temperature dependence of the resonance field H z for a 
given k gives the information on the leye]— a that domi- 
nates tunneling at a given temperature. EjE3 Eq. J72]) re- 
mains correct if the resonances for different n take place 
at different values of the field H z , provided that the field 
is being swept through all these resonances. If different 
n-resonances are well separated from each other, which is 
the case for Mni2 , then for each n-resonance one obtains 
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(73) [cf. Eq. (|7 



where Se n k is the level spacing in the well. One obtains 
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VII. LANDAU-ZENER RELAXATION IN 
SYSTEMS WITH DISTRIBUTED PARAMETERS 

As follows from the above, in M1112 crystals with dis- 
locations the tunnel splittings A„fe differ from one Mni2 
molecule to another. Thus the fraction of Mni2 molecules 
that stay in the metastable well after crossing the n- 
resonance, 
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is an average that must be computed with the help of the 
distribution function of Eq. ( |6l| ) 
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wher e f n k is given by Eq. (|76|). As was argued after 
Eq. (plj), the distribution of splittings is so broad that 
the natural scale to represent it is logarithmic. On the 
logarithmic scale, one can replace the exponential by a 
step function, e x => 8(1 — x). This gives 
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where A 2 k v is the solution of the equation 

vr/ AL (A 2 )/(25) = l 
for A 2 . With the help of Eq. @ we get 

'5s nk 



A 2 , 

nk.v 



00 for v > Wmax, where 

-max _ Se n k 



-Ae„ fc /T 



2tt 



(80) 



(81) 



(82) 



(83) 



The physical meaning of Eq. ( [79[) is the following. 
Since the tunnel splitting A is a high power of the trans- 
verse anisotropy E, there is very few Mni2 molecules 
with A n k ~ A. n k,v> whereas most of the molecules have 
A rlfc < A„k, v or A„ fc > A„k,v Eq- (© simply ignores 
the small group of Mni2 molecules with A n k ~ A n k }V - 
The quantity w™ fc ax is the maximal sweep rate v for which 
the Landau-Zener transition for the given n, k, and T is 
possible. For v > w™ fc ax one cannot find Mn 12 molecules 
that relax fast enough to satisfy Eq. ( |79|) , since the tran- 
sition probability has an upper bound on A which is de- 
scribed by Eq. (|76|). Eq. ( |S0f ) can be rewritten in terms of 
the distribution function for the logarithm of the trans- 
verse anisotropy, Eq. ([56]). It then becomes 



Rnk — 



dxf L (x), 



(84) 



where x n k. v — ln^nfc,« an< i E n k, v is the transverse 
anisotropy needed to create the value of splitting A^^. 
With the help of Eqs. © and @ we get 
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(85) 



One can see from Eq. ( p4| ) that the relaxation curves 
Rnk for the n, fc-resonance, when represented in terms 
of x n k,vi depend only on the distribution of transverse 
anisotropics in the crystal and scale onto the same 
curve. If the Gaussian approximation for the anisotropy- 
distribution function is used, see Eqs. (|53|) and (|57|), 
Rnk{xnk,v) can be expressed via the error function. The 
relaxation curves at T = are shown in Fig. ^| both in 
the scaling form, as functions of —x v , and in the nat- 
ural form, as a function of log 10 (l/z;) for even k. For 
odd k the relaxation curves (not shown) differ by a shift 
from the ajacent even-fc curves, see Fig. ^. The curves 
in Fig. ^|b have been plotted with the help of the simpli- 
fied Eq. (p0|). The deviation from the exact Eq. ( [79^ is 
rather small and it can be further reduced if one uses the 
improved formula 



R>, 



dxf L (x) 



C 



k.v 



-fL{Xnk,v), (86) 



where C — 0.577216 is the Euler constant. Note that 
the 1 / correction above has a non-scaling form. This 
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correction is equivalent to a small shift of the relaxation 
curves in Fig. [)] to the left: 



X n k =^ X n k — X n k C/(2f n J ). 



(87) 



As we have seen in Sec. distributions of tunnel split- 
tings for any n and k are functions of the distribution of 
transverse anisotropics, which is the basic characteristic 
of the Mni2 crystal. Field-sweeping measurements in the 
subkelvin temperature range, where transitions via ex- 
cited levels, n > 0, are negligible, provide the means of 
extracting that distribution. For odd fc-resonances, one 
should use Eq. (^34j) and fit the constant c oc id- 

At higher temperatures, tunneling transitions via ex- 
cited states n > become activated. The resonances 
with higher n occur at lower values of the longitudinal 
field H z because of the small term — AS^ in the spin 
Hamiltonian. In an ideal Mni2 crystal, field-sweeping 
experiments could provide the information on the level 
n*(T) which maximizes the function /„& and thus dom- 
inates the escape process at temperature T. Here the 
field sweep should not be too slow, otherwise crossing 
the resonances with higher n will cause depletion of the 
metastable well and resonances with lower n will not be 
seen. In real Mni2 crystals with dislocations, the situ- 
ation changes drastically because of the broad distribu- 
tion of tunneling rates. The depletion becomes unavoid- 
able because for any v < w™ fc ax most of Mni 2 molecules 
in the crystal belong to one of two groups. Molecules 
of the first group cross the barrier with a probability 



close to one, while molecules of the second group stay in 
the metastable well with a probability close to one. The 
whole process strongly depends on the history of the sam- 
ple, in particular, on the direction of the field sweep. If 
the field is increasing, then at the n, fc-th resonance the 
Mni2 molecules which cross the barrier satisfy 



< A, 



but 



< A 



(88) 



where A n k,v is given by Eq. ( |82[ ). The molecules with 



< A n+ x,jfc cannot cross the barrier at the nth 
resonance because they have already left the metastable 
well at the (n + l)-th resonance. The inequalities in Eq. 
(p8|) determine the range of the logarithm of the trans- 
verse anisotropy x = In E 



(89) 



for molecules that are going to escape at the nth reso- 
nance. The fraction of these molecules, which determines 
the magnetization step at the resonance, is given by 



dxf L (x). 



(90) 



If H z is increasing through a value at which the barrier 
disappears, all Mni2 molecules will escape, J2 n k ^nk,v — 
1. With increasing temperature or lowering the sweep 
rate, the depletion effects become stronger and stronger, 
which leads to the disappearance of tunneling resonances 
at certain n and k. Mathematically it manifests itself in 
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without hysteresis steps at n = resonances. 



the disappearance of the x- interval in Eq. (|89|), which 
leads to P n k,v = 0. 

To see which resonances are active and which are not, 
it is convenient to plot x n k, v of Eq. (|8^) versus n for 
different k, Figs. and [ll]. For dH z /dt > 0, the reso- 
nances are crossed in the order from small k to large k 
and, within a given fc-resonance, from the right to the 
left. If the value of x n k,v decreases, the x-interval in Eq. 
( |S9| ) does exist and the resonance is active. In the op- 
posite case the resonance vanishes. This happens for the 
resonances with high k in Figs. |l^ and [ll]. Note that 
the maximal value of n for each /c-resonance, n max , is the 
greatest n satisfying v < v™*. in the limit T — ► 
the maximal velocity u max of Eq. ( |S3| ) goes to zero for all 
n > 0. As a result, on lowering T, the branches of x n k, v 
in Figs. [To] and [ll] become shorter and shorter and even- 
tually reduce to the set of n = points along the vertical 
axis. This is in accordance with the obvious fact that at 
T = only transitions from the metastable ground state, 
ii = 0, are possible. With increasing temperature, the 
branches of move down for n > 0, and the ground-state 
resonances disappear one after the other. The same oc- 
curs if the sweep rate v decreases. The dependence of 
x n k,v on v, however, is only logarithmic [see Eq. ([S5l)1, so 
that the disappearance of the ground-state resonances at 
temperatures well below IK requires unrealistically small 
values of v. 

Note that the slope of the branches x n k,v in Figs. |l0| 
and [ll] changes for a certain combination of parameters k, 
T, or v for which the dependence x n k, v on n is nearly flat. 
One can estimate where it happens from the condition 
xok,v = x\k,v For the temperature T£(y) corresponding 
to this transition we get 



D(2S-\-k) 
\nQ v 



(91) 



This dependence is shown in Fig. |l2| for two different 
sweep rates corresponding to the boundaries of a typi- 



cal experimental window. In Mni2 crystals with disloca- 
tions, temperature T£(v) plays the role similar to that of 
temperature T q for an ideal Mni 2 crystal [cf. Eq. (6.1) 
of Ref . |l3| . The latter is the boundary between the pure 
ground-state tunneling and the thermally assisted tun- 
neling. It can be obtained by equating AL/fc fc and 
{A 2 lk /Se lk )e- A£lk / T . With the help of Eq. f|) one ob- 
tains 



T, 



00, k 



D{2S-l-k) 
hiQs 



Qi 



(92) 



In a crystal with dislocations the anisotropy E is dis- 
tributed and its role is played by the sweep rate v that 
"chooses" from the distribution of E the matching value 
of E satisfying Eq. ( |Sl| ) with n = 0. Substitution of the 
solution of that equation, Eq. (|85|), into Eq. (|9|), gives 
Eq. ©. 

The temperature Tqq^ is close to the quantum-classical 
transition temperature T (see Refs. [l3] |l4|) if the tran- 
sition is first order. Since for the model with transverse 
anisotropy the boundary between the first- and second- 
order transitions is Eq = D/3 (Ref. [l6]), which corre- 
sponds to Eq = 1/6, the values E <C 1 in the main part 
of the Mni2 crystal fall into the range of first-order transi- 
tion. Accordingly, the change between the two regimes in 
Figs. and [fl] is abrupt: The minimum of x n k v at n — 



is replaced by the minimum at n 



atT«T fe *({i) 



To study second-order transitions, one should use much 
greater sweep rates v which will shift the curves in Figs. 
[l0| and [ll] to higher values of transverse anisotropies E. 
This alone, however, does not guarantee that P n k,v of Eq. 
( pp| ) is large enough for the effect to be observed. For re- 
alistic concentrations of dislocations the values of /j, in 
that range of E are very small, see Fig. ||. 



VIII. RELAXATION IN THE THERMALLY 
ACTIVATED REGIME 

As was commented below Eq. (|76|), at temperatures 
above the quantum-classical transition, T > T the spins 
of Mni2 molecules escape over the top of the barrier 
through n = n top , since this value of n maximizes the 
function It is interesting to note, however, that 

Eq. ( |76| ) describes Landau-Zener transitions in the field- 
sweep setup, thus Eqs. (|74|) and ( f77| ) are only valid if it 
is the resonant tunneling that is dominating the escape. 
If the spins escape over the top of the barrier, it is no 
longer important whether the levels are in resonance or 
not, so that the steps in the dynamic hysteresis curves 
due to the Landau-Zener transitions should be washed 
out. 

Resonant transitions can still be detected at T > T 
due to the reduction of the effective energy barrier in 
measurements of the time relaxation and of the linear 
dynamic susceptibility because tunneling at resonance 
occurs via the level pair for which A n k is comparable 
with the level width T nn > (Ref. |l3|). These levels are 
lower than those at the top of the barrier that satisfy 



i i 



, L 







H_ = 












c = 10" 2 , 10~\ 10~ 4 


10 4 


\ l f 3 




Exact 








Gaussian 















relaxation off resonance in the activation regime using 
the classical model, to demonstrate that the relaxation 
is non-exponential due to dislocations. 

For the classical model with transverse anisotropy E -C 
D, the energy barrier is given by 

U = U(E) S U - AU, 

U =DS 2 (l-h z ) 2 , AU= \E\S 2 {1- h 2 z ), (93) 
Ideal crystal 

where h z — H Z /[2S(D — \E\)]. If transverse anisotropics 
are distributed the magnetization relaxation curve in the 
Arrhcnius regime is given by 
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FIG. 13. Distribution of the transition temperatures Too 
of Eq. (^2|) in Mni2 crystals with dislocations. Dashed lines 
correspond to the Gaussian approximation of Eq. (|57|) 



A n k ~ Se n k- However, in the field-sweep setup the main 
contribution to the escape rate comes from the levels at 
the top of the barrier, which makes the escape process 
non-resonant. A seeming paradox is that resonant tun- 
neling transitions are observed in the dynamic hystere- 
sis experiments in the activation regime. □ This paradox 
can be resolved if one takes into account the distribu- 
tion of To in a non-ideal Mni2 crystal. Then some of 
Mni2 molecules have To < T, while others have To > T, 
making tunneling resonances weaker but still present as 
long as there are molecules with To > T. Since for real- 
istic concentrations of dislocations the typical values of 
the transverse anisotropy satisfy E D, the quantum- 
classical transition is first order and a good estimation 
for To is the temperature Too given by Eq. (f)2]). The 
distribution of T o can be obtained from the distribution 
of \nE : see Eqs. ( j56| ) and (|57|). This distribution has an 
appreciable width, especially for high concentration of 
dislocations, see Fig. In experiments, quantum steps 
in the hysteresis have been observed at temperatures as 
high as 2.8K. This may be either an indication that the 
fine tuning of the theory (see below) is needed or an indi- 
cation that Mni2 crystals, or even molecules themselves, 
contain some stronger defects than studied in this paper. 
Note that for E = En = D/3 (which is the boundary 
between the first- and second-order transition) one has 
T = SD/ttk 2K (Ref. ||). 

Since the first-order transition is sharp, there are two 
groups of Mni2 molecules: Molecules of the first group 
escape via thermal activation over the top of the bar- 
rier, while molecules of the second group escape via the 
ground-state tunneling. Theoretical analysis of this situ- 
ation is more cumbersome since the basic Eqs. ( f74| ) and 
((77]) are valid only for a fraction of Mni2 molecules in the 
crystal. It is thus better to step back to Eq. ( |68| ) that 
yields a greater slope of the dynamic hysteresis curve at 
resonances and smaller but nonzero slope off resonances. 
We do not attempt to solve this problem in this arti- 
cle. Instead, we will consider for simplicity the time 



where 



/>oc 

R(t)=2 dEf E (E)exp[-r(E)t] 
Jo 



(94) 



(95) 



The transverse- anisotropy distribution fE(E) has been 
calculated in Sec. IV for a random array of linear dislo- 



cations. Using the results of that section, one can rewrite 
R(t) in the scaling form 



R(t) = 2 daf a (a) exp[-te po 



where t = t/r, r" 1 = Toe~ Uo/T , and 



P 



1 + h z 2U E C 



1 



T 



(96) 



(97) 



If p — > in Eq. ( pq ) , then one returns to the simple expo- 
nential relaxation, R(t) — exp(— t). In the case of p 3> 1 
one has R(t) = 2 f^" daf a (a), where a t = (1/p) ln(l/i), 
i.e., relaxation is logarithmically stretched. For Mni2 in 
the kelvin range, U /T is large but the value of E c is 
small for a realistic concentration of dislocations c, so 
that p is typically of order unity. Using estimates from 
the final part of Sec. IV, at H z = and T = IK one ob- 
tains p = 1.94 for c = 10~ 2 , p = 0.61 for c = 1CU 3 , and 
p = 0.194 for c = 10~ 4 . Thus Eq. @ does not simplify 
and it should be computed numerically with the use of 
Eq. (^7]). Note that for large crystals f a can be approx- 
imated by the Gaussian of Eq. ([49|), which implies that 
the actual parameter of the problem is p = 2\f~Lp w 6.5p. 

The short-time behavior of R(i) is singular due to 
Mni2 molecules [see Eq. (|52|)1 which are close to dislo- 
cations. These molecules correspond to the asymptote 
of the anisotropy-distribution function, f a = 1/a 3 , for 
a > «i, which yields 



R(t) = 1 



In(lA) 



(98) 



The singularity is rather weak for large crystals since the 
fraction n\ of molecules which relax according to Eq. d9§| ) 
is small. One can find the time t\ at which this stage of 
relaxation is completed from the condition R{t\) = 1 — n\ 
that gives 
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FIG. 14. Magnetization relaxation curves for Mni2 with 
dislocations in the activated regime. The fitted values of ( 
decrease on lowering temperature and increasing c. 



ti = exp[-p/y/nZ\ = exp 



4L ln(4L/-v/7r) 



(99) 



It is exponentially small for p 3> 1 and/or for large crys- 
tals. If [see comment after Eq. ( p2|) one replaces a% by 
three 3ai, the exponent of Eq. (|99) will increase by a 
factor of three. This means that the actual value of t\ is 
much smaller and the best way to find it is to c omp are Eq. 
( |98| ) with the exact relaxation curve of Eq. (M). Note, 
however, that in natural units this time is t\ = rti, where 
r is the relaxation time of an ideal Mn 12 crystal which is 
exponentially long for p> 1. 

Numerically computed relaxation curves R(t) for H z = 
0, T = 2K, and different concentrations of dislocations c 
are shown in Fig. [u| One can see that deviations of R(t) 
from a pure exponential are quite pronounced. For not 
too large concentrations c one can fit R(t) by a stretched 
exponential 



R(t) = exp(-at c ). 



(100) 



Although this dependence does not follow from any the- 
ory, the fits are surprisingly good. It would be i ntere sting 
to see if the fit of the relaxation curve with Eq. (100) pro- 
vides the actual concentration of dislocations in a Mni2 
crystal found by other methods. For this purpose, we 
listed the fitted values of a and £, together with the pa- 
rameter p of Eq. (p7[), in Table | for Mni2 crystals of the 
typical size 0.5 x 0.5 mm 2 at T = 2K and H x — 0. For 
other temperatures and magnetic fields, one can interpo- 
late a and £ on p using Table |. Note that in our picture 
the exp(— y/t) relaxation observed in some experiments is 
one of many possibilities corresponding to different fields, 
temperatures, and concentrations. 



TABLE I. Parameter p of Eq. @ and fitting parameters 
a and £ for 0.5 x 0.5 mm 2 Mn u crystals at T = 2K, H z = 0, 
and different concentrations of dislocations c. 



c 

3 x 10 -4 
10 -3 

3 x 10~ 3 
10 -2 



P 



0.097 
0.168 
0.306 
0.530 
0.97 



1.38 
1.68 
2.19 
2.68 
3.13 



0.97 
0.91 
0.77 
0.59 
0.40 



into account imperfections of the crystal structure and is 
based upon the generalization of the Landau-Zener ef- 
fect for incoherent tunneling from excited energy levels. 
All experimental features of the low-temperature mag- 
netic behavior of Mn 12 crystals find natural explanation 
within this theoretical framework. 

Linear dislocations at plausible concentrations provide 
the transverse anisotropy which is the main source of tun- 
neling in Mni2. Local rotations of the easy axis due to 
dislocations result in a transverse magnetic field for any 
external field applied along the c-axis of the crystal. For 
odd resonances, the tunneling matrix element contains 
the transverse field only in the first order of the pertur- 
bation theory, while the tranverse anisotropy makes the 
principal contribution to the transition amplitude. This 
explains the presence of odd tunneling resonances and 
their strength relative to that of even resonances. 

Transverse anisotropies and fields are distributed in the 
crystal. One consequence of that distribution is that the 
temperature of the crossover between quantum and ther- 
mal behavior becomes distributed within an appreciable 
range. This may be in part responsible for the experi- 
mental fact that quantum steps in the hysteresis of Mni2 
are well pronounced in the thermally activated regime. 

Crystal defects produce a broad distribution of tun- 
nel splittings that can be extracted from the dependence 
of the magnetic relaxation on the field-sweep rate. The 
theory predicts that relaxation curves for different tutu 
neling resonances can be scaled onto one master curve.EEl 
The first derivative of this curve equals the distribution 
function of transverse anisotropies in the crystal. 

Due to that distribution, quantum steps in the magne- 
tization curve should appear and disappear in a peculiar 
manner. At a given temperature and field-sweep rate, our 
theory predicts which steps must be abseniO regardless 
of the distribution function. 

Another consequence of the distribution is that the 
magnetic relaxation in the thermally activated regime fol- 
lows the stretched-exponential law. The exponent in this 
law depends on the field, temperature and concentration 
of defects. In zero field at T = 2K it is between 0.97 and 
0.40 for concentrations ranging from 10 -4 to 10~ 2 per 
unit cell of the crystal. 



IX. CONCLUSIONS 
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